June 2006 
TTP06-19 



Kl —> 7r°e + e and Kl —> 7r°/i + /i : 
A binary star on the stage of flavor physics 

Federico Mescia, 1 Christopher Smith, 2 Stephanie Trine 3 



1 INFN, Laboratori Nazionali di Frascati, 1-00044 Frascati, Italy 
2 Institut fur Theoretische Physik, Universitat Bern, CH-3012 Bern, Switzerland 
3 Institut fur Theoretische Teilchenphysik, Universitat Karlsruhe, D-76128 Karlsruhe, Germany 



Abstract 

A systematic analysis of New Physics impacts on the rare decays Kl — ► tt°1 + 1~ is performed. 
Thanks to their different sensitivities to flavor-changing local effective interactions, these two modes 
could provide valuable information on the nature of the possible New Physics at play. In partic- 
ular, a combined measurement of both modes could disentangle scalar/pseudoscalar from vector 
or axial-vector contributions. For the latter, model-independent bounds are derived. Finally, the 
Kl — > 7r°/i + /x~ forward-backward CP-asymmetry is considered, and shown to give interesting 
complementary information. 



1 Introduction 



Rare K decays, directly sensitive to short-distance FCNC processes, offer an invaluable window 
into the physics at play at high-energy scales. Besides the two K — > ttuu golden modes, the decays 
Kl — > 7r°e + e~ and ifx — ► 7r°/i + /i~ also exhibit good sensitivities, thanks to the theoretical control 
achieved over their long-distance components [1-3]. Importantly, these modes are sensitive to 
different combinations of short-distance FCNC currents, and thus allow in principle to discriminate 
among possible New Physics scenarios. 

In this respect, the pair of Kl — > tt°£ + £~ decays is unique since, though their dynamics is similar, 
the very different lepton masses allow to probe helicity-suppressed effects in a particularly clean 
way. Only the Kl — > £ + £~ modes share this characteristic, but the dominance of the long-distance 
two-photon contribution unfortunately prevents from acceding to the short-distance physics with 
a good degree of precision [4]. On the contrary, the corresponding two-photon contribution to 
Kl — > tt°£ + £~ is under control. It represents only 30% of the total rate for the muonic mode, and 
is negligible for the electronic one [2,3]. 

The main purpose of the paper is to illustrate how this fact can be used to constrain or identify 
the nature of possible New Physics effects. More precisely, our goals are: 

1. To analyze the impacts arising from all possible AS = 1 four-fermion operators in a model- 
independent way, i.e. operators of the form (sFd) (£T£), T = 1, 75, 7**, 7^75, , and to 
show how combined measurements of Kl — ► ir°£ + £~ can disentangle them. An important 
distinction is made between helicity-suppressed operators, like the T = 1,75 ones arising 
for example in the MSSM at large tan/3 [5,6], and helicity-allowed operators like in SUSY 
without R parity [7] or from leptoquark interactions [8]. Also, the electromagnetic tensor 
operator sa^ u dF^ u will be briefly considered [9, 10]. Finally, constraints from Kl — * £ + £~ 
(for scalar /pseudoscalar operators) and K — > 'kvv (for helicity-allowed tensor /pseudotensor 
interactions) will be analyzed. 

2. To improve the control over the long-distance two-photon contribution. Indeed, this is needed 
to estimate interference effects with New Physics short-distance contributions. In addition, 
once achieved, the Kl — ► vr°/i + /x~ forward-backward (or lepton energy) CP-asymmetry [11,12] 
will be computed reliably for the first time, both in the Standard Model and beyond, and will 
prove to be an interesting complementary source of information on the New Physics at play. 

The paper is organized as follows. In Section 2, we present the ingredients needed to deal with 
the long-distance dominated contributions, and, in the spirit of Ref. [3], analyze the possible signals 
of New Physics in the vector or axial-vector operators. Then, in Section 3, we analyze all other 
operators, both in the helicity-suppressed and helicity-allowed cases. The corresponding analysis 
of Kl.s £ + £~ is in appendix. Finally, our results are summarized in the Conclusion. 



2 Kl —> tt £ + £ with standard short-distance operators 

The Kl — > ir°£ + £~ decays receive essentially three types of contributions, depicted in Fig^ 

A first class of effects, purely sensitive to short-distance physics, results from heavy particle 
FCNC loops (the W and Z bosons and the t and c quarks in the SM, see FigJ^i), and can be 
parametrized by a set of local effective operators. In the SM, the leading relevant effective Hamil- 
tonian induced by these effects reads [13]: 

<// = -%T^ [y 7V (s lfM d) (I-fi) + yr A (s^d) [hV)] + h.c. , (1) 
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Figure 1: (a) Short-distance penguin and box diagrams for the initial conditions of Vtajv i n the 
SM. (b) K° — K° mixing-induced contribution, with a long-distance dominated, CP-conserving 
effective K$ — ► 7r°7* vertex, (c) Long-distance two-photon induced CP-conserving contributions. 



with a = a (Mz) and X g = V* s V q d- Of course, in the presence of New Physics, other types of 
effective interactions could be produced. This will be the subject of Section 3. For now, we will 
assume that New Physics affects only the values of the coefficients y-j a,7Vj an d leave them as free 
parameters [3]. 

Beyond SM scenarios leading to such modifications of the vector and axial-vector couplings 
are numerous. Examples are the MSSM for moderate values of tan/3 (see e.g. Ref. [14], and 
references therein), for large tan /3 (from charged Higgs penguins, see e.g. Ref. [15]), or the Enhanced 
Electroweak Penguins of Refs. [16, 17]. Of course, in specific models, New Physics also affects 
operators with different flavor quantum numbers. We have opted here for a decoupled, model- 
independent analysis, considering thus only the operators relevant for Kl — ► tt £ + £~. 

Also, the four-quark operators Qi t ...fi have not been explicitly included in Eq.Q. This is 
because their impact on the direct CP-violating (DCPV) contribution to Kl — ► tt°£ + £~ , associated 
to the local effective Hamiltonian 7~L^ff m standard terminology, can be safely neglected in the 
SM [2,13]. New Physics cannot change this picture as new sources of CP-violation from Qi ... e are 
bounded from purely hadronic K decay observables. 

A second class of contributions, dominated by long-distance dynamics, is driven by the coupling 
of leptons to photons, via mixing (Fig^>) or via a two photon loop (FigJ^). Let us now 

analyze these in more detail. 

2.1 Long-distance dominated contributions 

The remarkable point with the contributions depicted in Fig. lb and lc is that they can be entirely 
determined from experimental data. Their estimation is thus not affected by possible New Physics 
effects. As they remain as an unavoidable background to the interesting short-distance contribu- 
tions, we briefly recall (and partly improve, in the case of two-photon amplitudes) the way they 
are dealt with. 

The indirect CP-violating contribution (ICPV, Fig^») originates from K° — K° mixing. 
The subsequent CP-conserving K$ — ► 7r°£ + £~ decay is dominated by the long-distance process 
Ks — ► ti" (7* - ► producing the lepton pair in a 1 state. The corresponding amplitude 
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reads: 

M l cPV = eM{K 1 ^n £ + r)=-e^^Ws(z)(P + Kr{u pW } , (2) 

where P,K,p and p' denote the momenta of the K,tt,£~ and £ + states, respectively, z = T 2 /M^o 
with T = p + p' and a em ~ 1/137. The Ws function has been analyzed in detail in Chiral 
Perturbation Theory (ChPT) in Ref. [1], and can be parametrized as follows: 

W s (z) = a s + b s z + W]F (z) + W$ K (z) . (3) 

The pion and kaon loops (W^ ,W{f K ) were found small and, to a good approximation, a single 
(real) counterterm dominates: Ws (z) ~ as- This counterterm can then be extracted from the 
experimental Ks — > 7r°e + e~ and Ks — ► 7r /U + /i~ branching fractions: \as\ = 1-20 ± 0.20 [18]. 
Eq.(J2J) is thus indeed entirely determined in terms of measured quantities (\as\, e). 

T/te two-photon CP-conserving contribution (Fig^), Kl — ► 7r° (7*7* — > £ + ^~), produces 
the lepton pair in either a phase-space suppressed tensor state 2 ++ or a helicity-suppressed scalar 
state ++ . The former is found negligible from experimental constraints on Kl — ► 7r° (77)2++ [2], 
while for the latter the amplitude reads: 

M, n = ^0n MK r j^ (z) {u p v p ,} , (4) 

with 77 = raijMx and |Grg| = 9.1 • 10~ 12 MeV~ 2 . It is dominated by the two-loop process 
Kl — * vr° (P + P~ — * 7*7* — ► £ + £~) 0++ with P = n,K, computed in ChPT, and closely related 
to K s —> it + it~ — * 7*7* — > With the parametrization .M (-Ki — * tt°P + P~) ~ af (z) for the 

momentum distribution entering the subprocess P + P~ — > 7*7* — > the two-loop form-factor 

^ (2) can be expressed as (r„- = M^/Mk) 

F e ( z )=al(z)l( r lA)-a?(z)l( r l±) , (5) 



z z \ z z 



with I(a,b) given in Refs. [3,19] (for practical purposes, a numerical representation is given in 
Appendix B). 

A reliable estimation of T {Kl — * tt°£ + £~^ can then be obtained from the measured Kl — > 
7r°77 rate [20] thanks to the stability of the ratio 

r {K L -» ir°£+£-)^ 
R ^n = riKL^n^)^ ' (6) 

with respect to changes in af (z) [3]. Note that some O (p 6 ) ChPT effects are thus included in the 
estimated V {Kl — > 7r°.£ + ^~)^, most notably those responsible for the large observed Kl — ► vr°77 

rate compared to the (p 4 ) ChPT prediction. 

However, theoretical control over i? 77 only is not sufficient to deal with New Physics interactions 
that produce the lepton pair in a ++ state (generating interference effects) or compute forward- 
backward asymmetries. For this, we need to control M.~ n too. To fill this gap, the key is to use the 
similarity of behavior of the Kl — ► 7T°77 and Kl — > ~k°£ + £~ spectra at the origin of the stability 
of Rrpy. First consider the fact that the Kl — ► 7r°77 normalized spectrum is rather well-described 
with the parametrizations 

O [p A ) ChPT : a\ (z) ChPT = z - r 2 , (7a) 
K L -> vrV+vr" Dalitz : of {z) Dalltz = -0.46 + 2.44z - 0.95z 2 , (7b) 
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and (z) = z — 1 — r 2 in both cases. For a\ (z) Dalitz , only the z-dependent part of the Kl — ► 
7r°7r + 7r _ effective vertex is kept. Now, to account for the rescaling of B [Kl — > 7r°£ + £~)^ by 

B (K L -> 7r°7 7 ) exp , we rescale of (^) ChPT -» 1.45af (^) c/lPT and of (z) Dato2 -» 1.36of (^) Da ; ite , 
such that the rate computed from Eq.@ is kept frozen. Clearly, the theoretical control on the 
resulting Kl - > ir°l + £~ differential rate is not as good as on the total rate, but is nevertheless 
satisfactory. In practice, the error inherent to the procedure can be probed by comparing the 
predictions obtained using either a\ (z) ChPT or af {z) Dalitz . 

Finally, as far as the sign of .A/f 77 is concerned, we checked that Eq.@ is consistent with 
the conventions used in the rest of the paper for hadronic matrix elements. Furthermore, under 
the reasonable assumption that the sign of G% as fixed by the factorization approximation is not 
changed by the non-perturbative evolution down to Mk (see for example Ref. [21]), our conventions 
correspond to Gg < 0. 

2.2 Vector and axial- vector short-distance contributions 

Let us now consider the DCPV piece induced by the vector and axial- vector operators of Eq.Q: 

-Mv,A = (n°t + r\ -H^j\K L ) = iG F a(n°\s^d\K ) {u pl ^ (lm(X t y 7V ) + lm(X t y 7A ) 75 )v pl } (8) 

(the sizeable c-quark contribution, known at NLO [13], is understood in y 7 v), with the matrix 
clement 

(vr (if) \irfd\K (P)> = -j= ((P + Kf ff*° (z) + (P - Kf f«°*° (z)) , (9) 

f«°*° (*) = ^ (ff*° (z) - /f"° (z)) . 
The slopes are extracted from K ^ decays [22,23], neglecting isospin breaking: 

/„,+ (z) = f«l«° (z) = / + . (0) , A = 0.18, A + = 0.32 , (10) 

in the pole parametrization. Accounting for isospin breaking in ir° — r\ mixing for the value at zero 
momentum transfer [24], one gets: 

/+ (0) = (1.0231) ff V+ (0) « 0.939 , (11) 

with the Leutwyler-Ross prediction /^ 0?r+ (0) = 0.961(8) [25], confirmed by lattice studies [26]. 
This quite precise value, together with the knowledge of the form-factor slopes, renders the theo- 
retical prediction of the vector and axial-vector contributions to Kl — ► ir°£ + £~ remarkably clean. 

Inserting Eq.@ in -Mv,A, the vector current is seen to produce the lepton pair in a vector state 
1 , while the axial-vector part produces it in an axial- vector 1 ++ or pseudoscalar ^ state. This 
latter component is helicity suppressed, and enters thus only the muon mode. Besides, the vector 
current and ICPV amplitudes produce the same final state, they thus interfere in the rate. Recent 
theoretical analyses point towards a constructive interference, i.e., as < [2,27]. 



Total rates: Altogether, the branching ratios are predicted to be 

&V,A = ( C dir ± C int \ a s\ + C l mix |a 5 | 2 + • 10" 12 , (12) 

C% r = (4.62 ± 0.24) {w 2 7V + w 2 7A ) , C* r = (1.09 ± 0.05) {w 2 7V + 2.32w 2 7A ) , 

Cf nt = (11.3 ± 0.3) w 7V , Cf nt = (2.63 ± 0.06) w 7V , 

C^ = 14.5 ±0.5, C^ = 3.36 ±0.20, 
C^raO, = 5.2 ± 1.6, 
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Figure 2: Behavior of B {Kl — ► 7r °A t+ A i_ ) against £> (iTj, — > 7r°e + e _ ), in units of 10 , as ?/7A,7V 
are rescaled by a common factor (dots) , or allowed to take arbitrary values (red sector) . The ellipses 
denote 25, 50, 75% confidence regions in the SM, assuming constructive DCPV - ICPV interference 
(destructive interference is around the —1 dot), or for Ref. [17], in which uta ~ —3.2 and y^y ps 0.9. 



with wtajv = I m (-^2/7A,7v) /ImAt. In the Standard Model, the coefficients y7A,7V are real [13]: 



y 7A (M w ) = -0.68 ± 0.03, y 7V (jm ps 1 GeV) = 0.73 ± 0.04 , 
and, with lm\ t = (1.407 ± 0.098) • 10~ 4 [28], one gets 

Btl~= 3.54^1 (l.56lgffl ' ID" 11 , B^" = 1.4118^ (0.9Big^) • 10" 



ii 



(13) 



(14) 



for constructive (destructive) interference. The present experimental bounds are one order of 
magnitude above these predictions: 



B e e ~ 

*^exp 



< 28 • 10 
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29], jB£T 



< 38 • 10 
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[30]. 



(15) 



The differential rate for the electronic mode is trivial (no CPC piece), while for the muonic mode 
it can be found in Ref. [3]. 

The Standard Model confidence region on the B e e —B^ ^ plane is shown in Fig|2J As dis- 
cussed in Ref. [3] , this plane is particularly well-suited to search for New Physics signal, and identify 
its specific nature. Indeed, the fact that helicity suppression is rather inefficient for the muonic 
mode introduces a genuine difference of sensitivity to the short-distance V,A currents, i.e. to y-jA^v-, 
for the two Kl — ► tt°£ + £~ modes. These two types of contributions can thus be disentangled by 
measuring both rates. Suffices to note that the coefficients approximately obey Cf/C? ps 0.23, 
i = dir,int,mix, which is simply the phase-space suppression, except for the enhanced yjA contri- 
bution to Cj- , which comes from the production of helicity-suppressed pseudoscalar h states. 
Without this, no matter the New Physics contributions to y7A,7V, the rates would always fall on a 
trivial straight line. Thanks to this effect, on the contrary, for arbitrary values of y7A,7V, the two 
modes can lie anywhere inside the red sector in Fig E] for the Cf at their central values. Accounting 
for theoretical errors at la, this area is mathematically expressed as 



0.1 • 10~ n + 0.24B^'£ <B£X < 0.6- lO" 11 + 0.585 



-ii 



V,A 



(16) 
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Figure 3: Left: Ap B (in %) as a function of z, in the SM. Right: Integrated A FB (in %) as a function 
of wjv, with wja fixed at its SM value. Red (blue) lines correspond to constructive (destructive) 
interference between ICPV and vector current contributions, while plain (dashed) lines correspond 
to a i ( z ) D aiitz K ( z )chPT)i respectively. 



For definiteness, we have also indicated the curve corresponding to a common rescaling of yjA and 
y-jy from New Physics, and drawn the confidence region for the enhanced electroweak penguins 
of Ref. [17], to illustrate the opposite situation in which yjA = —3.2 is strongly enhanced, while 
U7V = 0.9 stays roughly the same as in the SM. Finally, note that the extent of the confidence 
regions essentially reflects the uncertainty on as (whose effect is included in the bounds Ea.l|16|)l. 
and could be reduced by more precise measurements of the K$ — ► tt 1 + £~ branching fractions. 

Forward-backward asymmetry: The differential forward-backward asymmetry [11, 12] is de- 
fined by 

r»> d 2 T f° d 2 T 

- d y - / -r-r d y 



with 



A t f Js_ ±o d y dz J-v d y dz h7 . 



> i d y + i~i~ d y 

o dydz J_ yo dydz 
p . ( v — j/) i 

y= m 2 ' w = 2 AAr> (18) 

K 

and Pj = 1 — 4rj/z, f3 2 = A (l, r 2 , z) , A (a, b, c) = a 2 + b 2 + c 2 — 2 (ab + be + ac) . The variable y is 
related to the angle between the K and £~ momenta in the dilepton rest-frame (hence the name 
forward-backward), and also, by definition, to the energy difference Eg- — E%+ in the K rest-frame 
(one then speaks of lepton energy asymmetry). 

This observable requires CP-violation, and arises from the Re[.A/^ (Ai ++ + .M 2 ++)] interfer- 
ence term. Let us start by assuming that M2++ 1S negligible, as for the total rates. Then: 

A'fb (*) = §J^ Z ^ (s*W* s (z) T t (z)) + Im (X mv ) hn^ (z) , (19) 

with A\ 2 (z) some combinations of constants and form-factors. The electronic asymmetry A e FB (z) 
is negligible in this case since J- e (z) is helicity suppressed, while for the muon it is shown in FigEl 
in the case of the Standard Model. Note that the theoretical control gained over this quantity 
would be difficult to improve since it relies on the specific parametrization of 7W 77 . In addition, 
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the experimental sensitivity required to measure A FB (z) is unlikely to be achieved soon. We will 
therefore not consider A FB (z) anymore, but rather concentrate on the integrated asymmetry: 

(l-M 2 / ryo d 2 T rO d 2 T \ 

dz / - — —dy — / - — —dy 

e _ Mr] \Jo dydz J- yo dydz J _ jy [E t - > E e + ) — N (E e - < E i+ ) 

FB- (1 _ Pir) a / ryo d 2 T ,o d 2 T x~ " AT (E t - > E e + ) + N (E t - < E e+ ) ' 1 } 

/ dz / tttt^ + / -r-r*/ 

y 4r | \Jo ay" 21 J- yo dydz J 

Compared to Ap B (z), it is more stable: 

A FB = (1.3 (l)w 7V ± 1.7(2) \a s \)-10- 12 / Bf£f . (21) 

In the Standard Model, ^4^ B = (20 ± 4) % for constructive and (— 12 ± 4) % for destructive inter- 
ference, with the error coming from varying a\ (z) between O (p 4 ) ChPT and Dalitz accounting for 
±2%. 

For general axial currents, it is clear that Ap B decreases when wia increases since it does not 
contribute to the interference Eq. QliJJI . For w-jy, the interference is linear while the rate is quadratic, 
and thus |^4^| reaches a maximum, around 23% for at its SM value, before decreasing again, 
see FigEl The absolute maximum for A FB is around ±25% for wia ~ and Wjv ~ if or =t4, 
depending on the direct - indirect CPV interference sign. 

In the SM, even if Ap B is polluted by the theoretical error on the two-photon amplitude, its 
measurement could fix the sign of 05. As can be seen in FigOU this remains true if New Physics is 
found from the measurements of Kl — ► ir°£ + £~ total rates with \wjv\ ^ 2. 

Let us now consider the interference term RefA'i^ M. 2 ++}- The amplitude Ai 2 ++ is discussed 

in detail in Ref. [2]. It leads to a helicity-allowed, but phase-space suppressed contribution to A e FB , 
hence contributes mostly for the electronic mode. Unfortunately, the theoretical control on the 
amplitude is not good as it depends on unknown phenomenological parameters. Though 
sufficient for deriving, from Kl — > 7r°77, a tight upper bound on the 2 ++ contribution to C^ 7 

in Eq.(|12|) [2], A e FB remains largely unconstrained. Indeed, both Ai* and Ai 2 ++ contribute 

mostly at low z, and their sizeable interference can generate A e FB anywhere between 0% and about 
±60% [11], depending on the phenomenological parameters of Ref. [2]. 

Concerning A„ B , the situation is better. The impact of A4 2 ++ corresponds to an additional 
±3% uncertainty, and therefore does not affect the potential of A FB in determining the sign of as- 



3 Kl — ► tt £ + £ with generic new physics operators 

In addition to the modification of vector and axial-vector couplings considered in the previous 
section, new four-fermion effective interactions could be generated by the integration of New Physics 
heavy degrees of freedom. The effective Hamiltonian comprising all the possible dimension-six 
semi-leptonic four-fermion structures [31] relevant for Kl — > 7r°£ + £~ (as well as quark bilinear 
electromagnetic couplings) reads 

n eff = n eff + n eff + Kff + H eff > ( 22 ) 
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with T~t e jf given in Eq.Q and 

oyP,s _ G F M w m s m e 
Tj f _ GpMyy m s m e 




[y P (id) (£ 75 £) + y s (Id) (#)] + fc.c. 



[ifr (S^d) (At* 4 "*) + y f (so^d) {ia^i)] + h.c. 



(23) 



(24) 



[y± (l<v (1 ± 75) d) F^] + h.c. . 



(25) 



A low scale (fi < m c ) is understood for the evaluation of the Wilson coefficients, quark masses 
and matrix elements of the operators in the above equation. Dimension-eight operators, containing 
two powers of the external momenta, are not considered as they are very small in the SM and are 
expected to remain so in the presence of New Physics (see discussion in [32]). 

Our goal is to analyze the impact of these new operators on the Kl - » tt®£ + £~ branching 
fractions and asymmetries in a model-independent way. Still, some comments on specific scenarios 
behind the various operators are in order: 

1. For the scalar and pseudoscalar operators Qs = (Id) (££) and Qp = (sd) (^75^), we have 
explicitly included the m s mi/M'y V helicity suppression factor to give a realistic description 
of models where these operators are generated from an extended Higgs sector. For example, 
large yps can arise in the MSSM with large tan/3 (see e.g. [5,6]) and sizeable trilinear soft- 
breaking couplings. This kind of scenarios has been analyzed in many works, but usually 
focuses on other decay modes. See e.g. Refs. [32, 33] for a MSSM analysis with emphasis on 
the K, B -» £+£~ decays. 

2. The tensor and pseudotensor operators Qt = {sa^ u d) (la^tj and Qf = (so^d) (ia^^i) = 
i E \ivpo (sa^d) (iupajbi) , to our knowledge, have not been included in studies of Kl — * 7r°£ + i~ 
so far. These modes are however the most promising source of information on Q T f, since these 
cannot contribute to K — * £ + l~ . Though they do not arise in the SM, they do in the MSSM 
but, in addition to being helicity suppressed, they are usually suppressed by loop factors [32]. 
Similar operators have been considered in AF = 2 processes (see e.g. Refs. [34,35]) and for 
the B — > X£ + £~ rate and asymmetries (see e.g. Ref. [36]). 

3. For completeness, we have included the dimension-five electromagnetic tensor operators = 
(scTfu/ (1 ± 75) d) F^ u . These were considered for example in Refs. [9,10]. Since the 07^75 part 
does not contribute to Kl — ► tt°£ + £~, only y y = y+ + y~ is accessible here. Note that in 
principle these operators already arise in the Standard Model, however they are too small 
to affect Kl — ► ir°£ + £~. In the MSSM, they are correlated with the chromomagnetic tensor 
operators and thus strongly constrained by other observables [9,10]. 

4. In the last section, we consider the general framework in which neither Qs,p nor Q T f are 
helicity-suppressed, i.e. we remove the m s mf/M^, factors in Eas. (|23l24|) . A large class of 
models with such helicity-allowed FCNC operators are theories with leptoquark interactions 
(for a review, see [8]), among which specific GUT models. Alternatively, SUSY without 
R-parity can also induce helicity-allowed Qs,p interactions through tree-level sneutrino ex- 
changes (see e.g. [7]). 

The distinction between helicity-suppressed and helicity-allowed scenarios is necessary as the 
corresponding signatures, i.e. impacts on Kl — ► 7r°e + e _ and Kl — ► 7r°/x + ^~, will obviously be 
very different. Let us now analyze these impacts systematically. 
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3.1 Scalar and pseudoscalar operators 

The relevant matrix element reads 

(Asd\K°) = - ^ K ~ M \ fo(z) (26) 
V2 {m s - m d ) 

in the sign convention of Eq.Q. This matrix element is enhanced compared to its vector coun- 
terpart due to the large value of the quark condensate (i.e., the large ratio of meson masses over 
quark masses). 

The scalar (pseudoscalar) operator produces the lepton pair in a CP-even ++ (CP-odd h ) 
state, therefore it is the real (imaginary) part of its Wilson coefficient that contributes to Kl — > 

(~<2 i\/[2 

•MP = « - Ml) l lmy P f (z) {u pl5 v p ,} , (27) 

^ = [M 2 R _ Reys /o {z) { - pVpi} (28) 

7r m w 

To reach these expressions, has been neglected against m s in Ea. H26() . The pseudoscalar current 
interferes with the helicity-suppressed pseudoscalar part of the axial- vector current, while the scalar 
current interferes with the helicity-suppressed two-photon ++ contribution. Since in addition Q$ 
and Qp are themselves helicity-suppressed, only the muon mode can be affected by ys,p- 

For the pseudoscalar operator, including also the V and A contributions, the differential rate 
reads: 

^ ViA&P _ G 2 F M\a 2 (Im At)2 ^ / Pl A+ {f+ {z)f + ^ (/q {z))2 \ j (29) 



dz 647r 3 n \ 6 

2 ' u z 



3 - ftp r 2 
A + = W 7A$ + W 7V — TT^i ^0 = — (w 7 A (l - rl) + zpplmyp)' 



with the prefactor 

1 1 Ml-M 2 

PP = —- = K 71 bs 0.18 . 30 

Im X t 2tt sm 2 6 w M% v 1 

Imyp ~ O(10) is thus required to get O (1) effects. Numerically, the contributions to the total 
rates, to be added to Eq.(|12jl. are 

Bf, +C ~ = (l.9w 7A Imy P + 0.038 (Imyp) 2 ) • 10~ 17 , (31a) 

B£ +fl ~ = (o.26w 7A Imyp + 0.0085 (Imyp) 2 ) • 10~ 12 , (31b) 

showing the very strong helicity suppression at play for the electron mode. 

For the scalar operator, from Eqs.® and l|28j). one immediately gets for the total ++ 
contribution: 

^r 77 fe g _ G\M\a\ m 3 rf 2 

— dz — — 5127r 7 — PnPe— \-r£{z) - ps&eyszfo {z)\ , (32) 



with the suppression factor 



V2tt aG F K-Ml _ 

PS = ^T79 ~ 0-13 • (33) 

sin 2 ^a2 m G 8 M 2 , 
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Enhancement 


Enhancement 


Maximal 


Bound from 


Experimental 




of 50% 


of 100% 


suppression 


K L -> 


bound (Eq.JTJJ)) 


Imy P 


-20, 40 


-30,50 


7% for Im yp « 10 


|Imyp| < 8 


|Imy P | <220 


Rey s 


-45,35 


-65,55 


1% for Reys » -5 


|Re2/ 5 | <50 


|Reys| <300 



Table 1: Numerical analysis of scalar and pseudoscalar operator impacts on B^ ^ . 



Since (z) ~ O (1), the helicity suppression M^y 2 turns out to be nearly compensated. Performing 
the z integral, we find 

Bf e ~ = (1.5 (3) Re y s + 0.0039 (Re y s f) ■ 10~ 16 , (34a) 
£>£f +M ~ = (0.04 (1) Re y s + 0.0041 (Re y s f) • 10~ 12 . (34b) 

The error on the interference term is estimated by varying the distribution aj (z) between O (jp \ 
ChPT and Dalitz (giving respectively 0.034 and 0.053 for the muonic mode). 

Total rates: The Qs and Qp operators do not affect the electronic mode due to their strong 
helicity suppression. For the muonic mode, combining Eas, H31l34|) with Ea.([12|) and fixing wja,7V 
at their SM values Eq. (jl3fl . the impacts on the total rate are summarized in the first three columns 
of Tabled 

A well- motivated scenario in which Reys and Imyp can be large is the MSSM for large values 
of tan/3 [5,6]. In that context, the contributions of Qs,p are related to those of Q' p = (sj^d) (J^sl) 
and Q' s = {s^d){ii): ys,p = y'psi with y' p and y' s further correlated. The contributions of Q' PS 
to Kl — > were analyzed for example in Ref. [33], with the result that values of a few tens for 

y' s P are compatible with AS = 2 and -B-physics data 1 . 

Without restricting ourselves to the MSSM, we can investigate the constraints on ys,p derived 
from the experimental Kl — ► rate under the assumption that ys,p = y'p s approximately 

holds. This analysis is presented in Appendix A, leading to |Imy s | < 8 and |Rey p | < 35 for 
y' 7A = yjA = —0.68 (for the Q' 7A = (sj^j^d) (^7^75^) operator). Allowing for New Physics in the 
axial-vector current with the bound \y' 7A \ < 3 corresponds to |Reyp| < 70, a much larger range 
since these two contributions interfere in the rate. 

Our numerical analysis can be summarized drawing the allowed regions on the B e e -BP * 
plane (Fig@J. Compared to the region spanned for general values of yjAjv (FigEJ) turning on 
ys,p basically extends the vertical spread since only the muon mode is affected, and this mostly in 
the upward direction (i.e., enhancements). This is illustrated in Fig0]by the light-blue/dark-blue 
region, corresponding to | Reys| < 90 / | Imyp| < 35, respectively. Imposing further the constraints 
from Kl — ► H + ^~ under the assumption y7A,s,P — v'-japs §i ves the yellow region. Obviously, the 
sensitivity of Kl — > 7r°/x + /i _ to Qs,p is quite good. 

It should also be noted that Qs,p contribute mostly for large z. While, as explained in Ref. [3], 
introducing a cut off at z ~ 4r| can reduce the contribution of the two-photon amplitude with 
respect to the Qtajv ones, thereby reducing the theoretical error, such a procedure would also 
reduce the sensitivity to Imyp and Reys significantly, and may thus not be desirable. 

1 Assuming a degenerate SUSY spectrum ~ ~ Ma ~ Mg, the gluino contributions to ys,p scale as [33] 

ys,p ~ {Mir /M\) tan 3 /3(l + 0.01 tan sign ^)" 2 + 18 (*Sb) ( 5 ") 32 ) ■ 

With tan (3 ~ 50, fi > (as favored from the muon g — 2) and Ma in the range 300 — 500 GeV [6], one gets ys.p ~ O (15) 
with the single mass-insertion {&ll(rr))h ~ 10 -2 (compatible with ek constraints [34]) and ys,p ~ O (30) with the 
double one (5rr) 13 { 5 ll) 32 ~ W~ 3 (constrained by AM Sjd [6]). 
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Figure 4: Left: Impacts of scalar and pseudoscalar operators in the B e e — B^ ^ plane of FigEJ 
Light blue (dark blue) corresponds to arbitrary uta,tv together with Re ys < 90 ( Im yp < 35), resp., 
while the yellow region corresponds to yjv,7A,s,P arbitrary but compatible with B (Kl — > / u+ A i_ ) GXP 
(see text). The dashed, light-blue line indicates the lower extent of the corresponding region. Right: 
The asymmetry A FB as a function of wjy and Re ys, assuming constructive interference (as < 0), 
y-jA = —0.68 and Imyp = 0. 



Forward-backward asymmetry: Reys enters the numerator of A FB through the interference 
term Ml CP yMs only: 

A FB = (1.3 (1) w 7V ± 1.7(2) \a s \ =F 0.057 \a s \ Rey s ) ■ W' 12 / ^,s,P • ( 35 ) 

There is no interference between .Ms and Aiy because of their 90° relative phase. For this reason, 
A FB goes to zero when y-jy and/or Reys becomes large, and reaches its maximum for moderate 
values, as shown in Fig|l]for yjA = —0.68 and Imyp = 0. If these latter two values are enhanced, 
since they contribute only to £>vASP' ^fb decreases, i.e. the figure remains the same but the 
absolute size of A pB is reduced. Finally, the figure for destructive interference is readily obtained 
by performing a vertical axis reflexion (w^y — > —wjy) followed by an overall sign change for A FB . 

No matter the New Physics behind y7A,7V,s,P, \A FB \ is always smaller than 30%, i.e., not very 
far from its SM value Eq. ([21|) . Given the theoretical errors, A FB does not appear very promising to 
get a clear signal of New Physics. Nevertheless, as said before, it offers a very interesting possibility 
of constraining the relative signs of Keys, yxv an d as when considered in conjunction with the 
K L -> ir°£ + £- total rates. 



3.2 Tensor and pseudo-tensor operators 



Let us now turn to the tensor operators of Eqs. (j24j) and (|25[). The relevant matrix element assumes 
the form 

p/^ K v — P u 

(n° (K) \sa^d\K° (P)> = i B T (z) , (36) 

with (7T° | sa^ v 75<i|if ) obtained through 0^75 = i£^ vpu a p(T . The tensor form-factor was studied on 
the lattice [37], with the result Bt (z) ~ 1-2/+ (0) /(l — 0.29z) at fj, ~ m c in the MS scheme (an 
earlier order-of-magnitude estimate may be found in Ref. [38]). 
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The electromagnetic tensor operator produces the lepton pair in a 1 state and the tran- 
sition is CP-violating: 

-Memo = i—j=- — lrD.y 1 B T (z) {u p pv p '} . (37) 

Let us take At = A+, which is good enough for our purpose, and can be justified in a pole model 
through the fact that z « M TV /M^ with Mt,v the nearest vector and tensor resonances. The 
effect of can then be absorbed into the vector-current Wilson coefficient w-jy [9] 

uw ^ v ^ = vw + —-—. (38) 



The above redefinition is independent of the lepton flavor, hence the Kl — ► tt°£ + £ modes cannot 



disentangle possible New Physics effects arising from from those arising in the vector current 



electroweak penguins and boxes. Such New Physics effects were analyzed in Section 2, see FigEJ 
The tensor operator also induces a CP-violating contribution: 

■Mt = Im A t flr (*) r iPT Im y T U p Un JP - P ' ^~ P ^ iy \ , (39) 

where we have defined 

= 1 2 m * M K „ 1 (m \ 

PT lmX t7rS m 2 e w ~ 4 1 Uj 

for m s Ri 150 MeV. As for the magnetic operator, the f part can be absorbed into wjy, but now 
the me dependence introduces an effective breaking of fj, — e universality in the vector current: 

w 7V -> Wj V = Wjy + p T \m.y T rj BT< f^ • ( 41 ) 

The second term in Eq . (|39j) is also CP-violating because P-(p—p') = yM\ is CP-odd. It produces 
the lepton pair again in a 1 state and can thus interfere with both the ICPV and vector operator 
contributions, producing an extra contribution to the A + factor defined in Ea. (|29l) : 

( A \ 2 B t(z) 2 , ,2 z Pf . 2 a 2 B r(z) ( i a e m Im (eW s (z)) \ 
( a +)t = r e - - n2 KPT Im yr) -^ + r e f3 e - p T Im y T w 7V - —= — - . (42) 



In the muonic case, the overall /3jj factor for the first term, corresponding to an orbital angu- 
lar momentum of two between the muons (y counts as one unit of angular momentum), makes 
this contribution significantly phase-space suppressed compared to the one shifting wjv, Ea.(|41[). 
Numerically, 

B^ e ~ = ((1(T 5 + 0.08) (lmy T ) 2 + ((29 + 15) w 7V ± (36 + 18) \a s \) Imy T ) • 10~ 19 , (43a) 
B^ + v~ = ((0.25 + 0.02)(Imy T ) 2 + ((29 + 4)u> 7 i/±(36 + 5)|a s |)Imy T ) • 10~ 15 , (43b) 



where the first numbers in each parenthesis come from Ea.l|41|). the second from Ea.(|42|). and the 
± sign corresponds to as = =F \a>s\- 
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The pseudotensor operator produces the lepton pair in a 1 + state, and is thus CP-conserving 

Mf = Im (z) r e p T Re yf — j^- P - {u p -y 5 v p > } . (44) 

As none of the other CP-conserving contributions produces such a final state, it represents a distinct 
contribution to the rate: 

^7 = ( ImA ^^ 2 (PTRey f ) 2 (B T (z)fz. (45) 
Numerically, this contribution is phase-space suppressed and very small: 

£| +e ~ = 7.9 (Rey f ) 2 • 1(T 21 , (46a) 

B^^~ = 4.9 (Rey f ) 2 ■ KT 17 . (46b) 

Total rate and forward-backward asymmetry: Overall, the effects of the Q T T operators are 
smaller than those of Qp.s because of the smaller matrix elements and the phase-space suppression. 
In addition, y TT < ys,p hi realistic scenarios, and tensor operators seem beyond reach. The only 
exception is the effective breaking of fi — e universality in the vector current induced by Imyx, 
which slightly extends downwards the "V,A only" region of FigEJwhen |Imyr| > 25. 

For Ap B , Im y<r enters through interference with Im (z) (not with Re ys since they are out 
of phase by 90°): 

A FB = (1.3 (1) w 7V ± 1.7(2) \a s \ T 0.057|a s | Rey s + 0.033 (4) lmy T ) ■ 10" 12 / B^ p T r (47) 

Qf does not contribute directly to A FB as it leads to a real amplitude, out of phase from the 
Qta one by 90°. Therefore, assuming Iinyy < Reys, no impact can arise for A FB . Similarly, the 
electronic asymmetry is not affected, even taking into account interferences with the M2++ piece. 

3.3 Helicity-allowed effective interactions 

We now lift the constraint of helicity suppression in Eqs. (|23|) and (|24j) . symbolically as: 



G 2 F M^m s m e g% P f 
More precisely, all former expressions remain valid provided one makes the substitutions: 

where the normalization Aj = (Aj/100TeV) corresponds to typical lower bounds obtained from 
lepton-number violating processes (see e.g. [39]), and we have assumed gNp/g 2 ~ Plugging 
these expressions in Eas. ()31I34I43I46|) . we obtain: 

+e - (115 2.6 112 0.3w 7A 2.3 0.07^ 0.09 \a s \ 




I 29 , 4.5 , 60 , 22w 7A , 1.9 , 2.8w 7V , 3.4 a s , , ,, 

B s^f=[^ + AT + A J + ^ + AT + ^ ± ^ + ^\- 10 • ( 50b ) 

Without helicity-suppression, it is the electronic mode that is the most sensitive to these operators, 
simply because of the phase-space suppression in the muonic mode. Therefore, allowing for these 
interactions typically produces total rates in the lower part of the B e+e -B^ + ^ plane, a region 
which cannot be attained by the New Physics interactions studied up to now (see Figs El and @J). 
Let us further investigate what signals could be expected. 
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Scalar /pseudo scalar operators: As before, assuming that the Wilson coefficients for the scalar 
and pseudoscalar operators contributing to Kl — ► £ + £~ and Kl — ► ir°£ + £~ have approximately 
the same values, Agr pi ~ Apg, the constraints from Kl — > /i + /*t~ still leave the possibility of 
sizeable effects. However, now that helicity suppression is no longer effective, one gets a very tough 
constraint from Kl — ► e + e~ . From the measurement B(Kl — ► e + e~) exp = 9^4 • 1CP 12 [22], and 
since 

B (K L - e+e") SiP » 6.8 • 1(T 8 ( ^- + , (51) 

one gets Agi, Ap/ > 8. For such large values, the effect on both is of a few percents for 2/7^4 at 

its SM value, hence well beyond reach. Though Agr pi and Ap t g can be different in specific models, 
the large difference needed to get observable effects on B e+e would require a somewhat fine-tuned 
scenario. We can therefore reasonably rule out this possibility. 

Tensor /pseudotensor operators: To get a somehow realistic estimate, we bound them from 
K + — > -k + vu assuming that the operators {sa fJ,v d)(va tll ,i') and (sa^ ' d)(va liv 'y^u) are governed by 
the same scale factors (gNP I 'A TT ) 2 . Their contribution reads 

B (K+ - x+vv) Tti « 4.4 • 10- 12 (j^ + , (52) 

and comparing with Eq.((2UJ) shows that the sensitivities of K — ► 7r^P and if^ - ► vr°£ + £ _ to tensor 
interactions are similar. 

Let us assume that A^ ~ Af,, and given the SM prediction of (8.0 ± 1.1) • 10" 11 [40] and the 
current measurement of 14.7^g 3 g° • 10 -11 [41] (V,A currents do not interfere with tensor operators 
for massless neutrinos), one finds A T j, > 0.35. Interestingly, for such small values, the charged 
lepton modes are quite large, B [Kl — * 7r°e + e~) as 3 • 10 -10 and B (Kl — > 7r°/^ + /i _ ) 2 • 10~ 10 , 
i.e., around their current upper limits Eq.(|15j). The conclusion of this order-of-magnitude estimate 
is thus that there is still room for large effects from tensor operators. The corresponding allowed 
region in the B e+e -B^ + ^ plane is shown in FigOD 

The forward-backward asymmetry: The muonic asymmetry A FB can be enhanced only by 
Qs and Qt contributions, as explained in earlier sections, and this through interferences with 
M 

icpv an d A4 77 _ ++, respectively, see Eq.(jHJ). Tuning y^y and the helicity-allowed 1/A T term, a 
maximum of about 60% can be reached. Still, this requires large contributions from Qt with both 
K L -> -K°e + £- rates above 10" 10 . With A T > 1, this maximum falls to about 30%, i.e. close to the 
SM value Ea.ljSTjl. 

For the electronic mode, one could think that a small Ag could generate a significant asymmetry 
through the helicity-allowed interference with A-^^py. However, this is not the case because the 
impact of Qs on B e+e is then much more pronounced. Tensor interactions, for their parts, have 
a small impact on A e FB because interference effects with M 2 ++ are helicity-suppressed. 

4 Summary and Conclusion 

The Kl — > tt°£ + 1~ modes offer a unique opportunity to probe a large range of New Physics AS = 1 
effective operators, see Tabled They are therefore an essential tool in the investigation of flavor 
structures beyond the Standard Model. Our work was to analyze, in a model-independent way, the 
possible experimental signatures of these New Physics interactions. 
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Figure 5: Impact of helicity-allowed tensor and pseudotensor operators in the B e e - B^ ^ plane 
of FigH] The green area (up to the green dashed line) corresponds to arbitrary y-jAjv together 
with A T j, compatible with B (K + — > n + vi>), as explained in the text. 



In the presence of vector and axial-vector interactions only, the two rates are bounded in the 
ge+e- plane ( gee FigEJ, i.e., at la, 

0.1 • 1CT 11 + 0.24B£X~ < Btff < 0.6 • 10~ u + 0.58E£a~ . (53) 

Any signal outside this region is an indication of New Physics FCNC operators of different structures 
(baring the possibility of large [i — e universality breaking in the V,A currents). We have identified 
two possible mechanisms. 

The first is from helicity-suppressed (pseudo-)scalar operators, as arising in the MSSM at large 
tan (3. These enhance the muonic mode without affecting the electronic mode. Such interactions 
would be revealed by measuring the two rates above the V,A region. Also, comparing with Kl —> 
shows, model-independently, that the Kl — > 7r°/i + M~ is more sensitive (besides being cleaner) 
to these types of interactions (see the yellow region in FigEJ). 

The second possibility is from helicity-allowed (pseudo-)tensor interactions, which could arise 
for example from tree-level leptoquark interactions. Because of the phase-space suppression, it is 
now the electronic mode which is more affected. These interactions would manifest themselves 
in a signal significantly below the V,A region (see the green region in FigEJ). Even assuming the 
presence of similar contributions for the K — * ixvv modes, there is at present no severe constraint 
on these effects. 

On the other hand, we found that both helicity-suppressed (pseudo-)tensor interactions and 
helicity-allowed (pseudo-)scalar interactions should not lead to observable effects. For the former, 
this is so because of significant phase-space suppression, while the latter are already strongly 
bounded by the very rare Kl — > e + e~ decay. Also, concerning the electromagnetic tensor operator, 
(sa^vd) F^, measurements of the Kl — > tt°£ + £~ rates cannot disentangle it from vector operator 
contributions. 
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Short-distance 


CP-property 


Helicity- Helicity- 




operator 


& J PC (££) 


suppressed allowed 


ICPV (K^-K ) 




CPV (1— ) 


Ea. ma 


Two-photon 




CPC (0++ 


En.(tT2l 


Vector 


\ if* / \ i j 


CPV 


Eq.(d2j) 


Axial- vector 


\ 1 f* / \ 1 1 <J / 


CPV (0~+ 1++) 


Eq.® 


Pseudoscalar 


(sd) {£ l5 £) 


CPV (0~+) 


Eq.(EU) Eq.jSUJ) 


Scalar 


(sd) (U) 


CPC (0++) 


Eq.(EU) Eq.® 


Tensor 


(sa^ u d) (ia^e) 


CPV (1~) 


Eq-gSl) Eq.® 


Pseudotensor 


(sa^d) (£a^ l5 £) 


CPC (1+-) 


Eq.(@nj) Eq.jSni) 



Table 2: Summary of the contributions entering the Kl — ► 7r°£ + ^ _ rate, and references to the 
relevant formulas in the text. The CP-property indicates which of the real or imaginary part of 
the respective Wilson coefficient contributes. Interferences occur whenever the li pair is produced 
in the same state. 'Helicity-suppressed/allowed' refers to the last four operators. 

For the (integrated) forward-backward asymmetry Ap B , the first reliable estimate has been 
obtained. It is typically of a few tens of percents, and can give important information, comple- 
mentary to the total rates. In the SM, it could fix the sign of the interference between the vector 
operator and ICPV contributions. Similarly, beyond the SM, it can be used to discriminate among 
various solutions once both K L -> ir°£ + £~ rates are measured. On the other hand, the electronic 
asymmetry A e FB is found either completely dominated by its (unknown) SM value, or too small to 
be of any use to constrain either New Physics or as- 

We have not included differential rates or differential asymmetries in the present study (though 
they can be trivially computed from our analyses). New physics does affect these observables, but 
they require a higher experimental sensitivity, so total rates and integrated asymmetries are more 
promising. 

With the general expressions for both Kl — ► ir°£ + £~ rates and asymmetries computed in 
Sections 2 and 3, the way is now paved for more model-dependent analyses. In this context, 
the B e+e -B^ +lJ " plane remains as a particularly convenient phenomenological tool to display the 
correlations among operators specific to a given model. 

In conclusion, the Kl — ► tt®£ + £~ system, together with the neutrino modes K — ► -kuv, has a 
considerable potential for unveiling/constraining the nature of possible New Physics in AS = 1 
FCNC, therefore playing an important role in the quest for a better understanding of the quark 
flavor sector. 
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A Constraints on (pseudo-) scalar operators from Kl — ► ri + ri 

We consider the following effective Hamiltonian: 

We// = ^ A ^A (*vy&d) (*y"V) + y (s^d) (Hi) + y' s ( S " 75 d) {££)] + fcc. 

(54) 

In the SM, y' 7A = y-jA and y' p s are negligible. Our goal is to get an order of magnitude estimate of 
the bounds set on y' p s by the measured Kl — ► fj, + fi~ rate. 
Using the matrix element parametrizations 

M 2 

(0|s 7m7 5^ (P)) = iV^FKP^, (0\s l5 d\K°) = -iV2F K , (55) 

m s + md 

in the same conventions as Eq.Q, the decay amplitudes (^ + ^~| — T~leff\KL,s) an d the total rates 
can be written as 

M (K L , S -> £+r) = 2 , ^ M ^ M % (^i i5 + S£, 575 ) V } , (56) 



r (tf Lt s -> £+r) = - 4r? ^(1 - 4r|) + B^ s j , (57) 

with As,Bl (Al,Bs) the CP-conserving (CP-violating) pieces given by 
M 2 M 2 

A{ = — f i Im y' s , B l L = — f Re - (2vr sin 2 %) (Re ( Vta) + *cVc) + , (58) 

M w w 

M 2 M 2 
4 = T#R«l/s + 4ry. ^ = v#iIm^-(2^sin 2 ^)zIm(A t y^) . (59) 



The c-quark contribution is negligible for K$ — * /J- + fJ- , while for Kl — > /i + // , it has been computed 
recently to NNLO giving y c = (— 2.0 ± 0.3) • 10" 4 [42]. Indirect CPV contributions are understood, 
M(K L ,s^i + n iCPV = eM(K h2 ^ £+£-). 



The two-photon term A s is given in [19] in terms of the two- loop form- factor of Eq.(J2J): 



A *~n = ~ 2G 2 M 2 \ (1 " * 1 (r '' * = 2 - 10 • 10_4 ( " 2 - 821 + ?L216) • (60) 



For Kl, the situation is less clear as the dispersive part of A L ^ is notoriously difficult to evaluate. 
Anyway, following the analysis of Ref. [4] , one can get the conservative estimate 



A ^ = V ^G 4 FlM^M| : 77) {Xdlsp + iXabs) = ±L98 ' 10 ^ ((0 ' 71 ± °- 15 ± L0) " ^ ' (61) 



The sign of this contribution depends on the sign of the Kl — * 77 amplitude [4] , itself depending 
on the sign of an unknown low-energy constant (see [21]). 

To get an order-of-magnitude estimate of the coefficients, we allow for both signs in the Kl — > 
branching ratio 

B(K L -» = (6.7+ (0.08Im^) 2 + (0.10 Re y' P + l.ly' 7A - 0.2 ± 0.4+°|) 2 ) • 10" 9 , (62) 

and reflect only the error on Xdisp- Then, imposing the rate to be within 3cr of the experimental 
value B (K L -> /x + /i"f xp = (6.87 ± 0.12) • 10~ 9 [22] corresponds to |Im^| < 8 and |Rey^| < 35 for 
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the SM value y' 7A = —0.68. Relaxing this latter constraint to \y' 7A \ < 3 corresponds to |Rey P | < 70, 
much larger since the two interfere in the rate. 

For Ks — ► , the experimental bound B (K$ — ► /t/ + / u _ ) cxp < 3.2 • 10~ 7 is still very far from 

the predicted rate of about 4 • 10~ 12 in the SM. Taking |Rey^| , |Imyp| < 50 cannot enhance the 
branching ratio beyond B (Ks — > ^ > _ ) < 10~ 10 . 

Finally, for completeness, the longitudinal muon polarization in Kl — > is expressed as [43] 

2 A /l-4r2lm(fl£X) 

i^L = — o o • (63) 

(1-4^)K| 2 + |^| 2 

In the SM, this quantity is entirely driven by the indirect CP-violating contribution, proportional to 
e, and is thus rather small, \Pl\ ~ 2 x 10 -3 [19]. Including the scalar and pseudoscalar currents, we 
find that for |Reyp| , |Imyp| < 75, only a 10% deviation of Pl from its SM value can be generated. 
On the other hand, |Rey^| < 75 can enhance \Pl\ up to about 10 -2 , while for |Imy^| < 8, \Pl\ 
can be as large as 8%. It is indeed this latter parameter which is the most important since it 
does not require any e factor. This makes Pl particularly sensitive to the presence of new CP- 
violating sources in the scalar operator (as discussed e.g. in [44]). Unfortunately, a percent level 
measurement of Pl in the medium term is unlikely, and the Kl — > 7r°/i + /i~ total rate is more 
promising to get a signal or set limits on these New Physics interactions. 



B Numerical representation of the two-loop form factor 

In Refs. [3,19], the two-loop form-factor is expressed as a complicated three-dimensional integral. 
For practical purposes, the following numerical representations can be used instead: 

r l r l\-{ ELo 4 (°- 33 - z )' k/2 > < z < °- 315 . 

i Efc=o&fc(^-0-30) _fc/2 , 0.315 O<0.6, 



with 



r 2 ^ 5 



z 



^ci(l-z)- fc / 2 , 0^z<0.6, 



k=0 



k 


a k 


K 







10.65 - il7.24 


54.76 - i23.60 


-26.069 - i3.6914 


1 


-6.964 -H13.31 


-53.89 - i37.02 


92.05 + i0.9337 


2 


1.312 - i2.521 


6.832 + i24.74 


-122.03 + i8.3238 


3 


-0.6045 + i0.3762 


0.942 - i5.509 


76.611 - i8.2021 


4 


0.113 - i0.0381 


-0.265 + ^0.5508 


-23.325 + ^3.0351 


5 


-0.0069 + i0.00186 


0.0154 - i0.0206 


2.7774 - i0.4036 




k 


A* 
a k 


K 


c k 





-8.5486 - i0.9184 


-5.0369 - il.4840 


-9.6037 + i3.7637 


1 


9.3942 - i0.4070 


2.8237 + ^2.9981 


26.408 - €14.04 


2 


-3.6104 + i0.8098 


-1.2503 - i0.7132 


-28.934 + ^19.15 


3 


0.7102 - zO.1910 


0.3303 + i0.1108 


16.159 - €12.15 


4 


-0.0696 + t0.0181 


-0.0412 - iO.0100 


-4.5365 + i3.715 


5 


0.00269 - i0.000577 


0.00194 + i0.000389 


0.5087 - €0.4428 



Using this parametrization, one can reproduce the rates and asymmetries typically with an error 
of about one percent, which is more than sufficient given the size of the theoretical uncertainty on 
the af (z) distributions. 
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